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c j I Abstract 

We demonstrate how the Moutard transformation of two-dimensio- 
na l Schrodinger operators acts on the Faddeev eigenfunctions on the 
zero energy level and present some explicitly computed examples of 
such eigenfunctions for smooth fast decaying potentials of operators 
with non-trivial kernel and for deformed potentials which correspond 
to blowing up solutions of the Novikov-Veselov equation. 

>: 

^ ; 1 Introduction 

1^ ■ In [I] we used the Moutard transformation for constructing examples 

of two-dimensional Schrodinger operators with fast decaying smooth 
■ potentials and nontrivial kernels and of blowing up solutions of the 

| Novikov-Veselov equation. 

y— I ■ Here we show how the Moutard transformation acts on the Faddeev 

eigenfunctions corresponding to the zero energy level. In particular, 

• i"H ■ we construct some explicit examples which correspond to operators 

y\ \ constructed in [1]. It appears that such Faddeev eigenfunctions have 

very interesting analytical properties and dynamics under the Novikov- 
Veselov flow. 

The presented procedure also enables us to extend in an essentially 
new way a list of two-dimensional Schrodinger operators for which the 
Faddeev eigenfunctions are explicitly computed at least for one energy 
level. 
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For multi-dimensional Schrodinger operator 



L = —A + u 

the Faddeev eigenfunctions ip(x,k),k £ C n ,x £ R n , are formal solu- 
tions to the equation 

Lip — Eip, E=(k,k), 

meeting the condition 

V> = e l{k ' x) (1 + o(l)) as|x|^oo. (1) 

Here (x,y) = X)J=i - T i2/i i s the standard inner product in C™. 

These functions were defined in [2] and play the crucial role in the 
multi-dimensional inverse scattering problem [3J 0] ■ 

In the sequel we consider two-dimensional Schrodinger operators: 

/ d 2 d 2 \ 
L = -A + u= - I + g^2 J + u = ~ 4dd + M ' 

where in the right-hand side we have d = = | — i-§^j ,d = 

m = h(m+ i 7k)> (^1^2) = {x, y) £ R 2 , z = x + ly £ C. 

In this well as for higher dimensions, there are not many 

examples for which these eigenfunctions were computed explicitly. In 
fact, they are as follows: 

1) for the point potential 

u(x) = eS(x), 

the Faddeev eigenfunctions are computed for all energy levels in [5]; 

2) for the Grinevich-Zakharov potentials introduced by Grinevich 
[B] and Zakharov (these potentials are reflectionless on this energy level 
E and decay as |^| -2 ) explicit formulas for these Faddeev eigenfunc- 
tions on the energy level E may be extracted from computations in 
0). 

In [7] it is shown that for negative energy levels E the functions 
x(x, k) = e~^ k ' x '>tp(x, k) = exp (—±y/—E (Az — j)) ip(x, k), where fci = 

(A + j) , fc 2 = — (A — j), satisfy the generalized Cauchy- 
Riemann equation 

d\ X ' 

i.e. they are generalized analytical functions, and if T is nonsingu- 
lar, then this level is below the ground state. Moreover in the latter 
case the inverse scattering problem of reconstructing potentials (with- 
out assuming that these potentials are sufficiently small) is solved by 
methods of generalized analytic functions theory [9]. 
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In this article we consider the Faddeev eigenfunctions on the zero 
energy level E = 0. They meet the following condition 



as \z\ — > oo. 



(2) 

In fact, this is one branch of these functions and another branch is 
given by solutions to the equation (—A + u)ip = with asymptotics 
e^ z (l + o(l)) as \z\ — > oo. For non-zero energy levels the analogous 
asymptotic expansions were given in [10] and for the zero energy level it 
is derived in a similar way from the analytical properties of the Green- 
Faddeev function of the Laplace operator. The detailed analysis of the 
latter function, in particular, was given in [5]. 



2 The Moutard transformation 

Given uj satisfying the equation 

Luj = (-A + u)uj = 0, 
the Moutard transformation of L is defined as 

L -> L = - A + u - 2A log uj = - A - u + 2- 

uj- 

It is easy to check that 

• if <p satisfies the equation Lip = 0, then the function 8 defined up 
to term —,C = const, from the system 

{uj6) x = -uj 2 (t\ , (uj6) y =uj 2 (^) (3) 

\UJ/y \UJ I x 

satisfies L9 = 0. 

If the potential u depends only on x, then a certain reduction of 
the Moutard transformation gives the famous Darboux transformation 
of one-dimensional Schrodinger operators: 

^ dx 2 \ dx W ) \ dx W ) ~~ ^ ^ \ dx ^ ^dx 

This is the case when u> = f(x)e y/By and + u^j f = Ef. 

In pQ we used a double iteration of the Moutard transformation 
starting from the Laplace operator Lo = —A. Given two harmonic 
functions uj± and 0J2, AcJi = Auj2 — 0, we construct a pair of Moutard 
transformations corresponding to uj\ and UJ2- 

Ly = -A + u x 4^- L = -A ^> L 2 = -A + u 2 . 
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Both potentials Ui — — 2Alogwi and U2 — — 2Alogu;2 have singular- 
ities at zeroes of functions u)\ and u>2, respectively, and we have to 
iterate the Moutard transformation to achieve smooth (i.e. bounded 
and diffcrcntiable) potentials. 




Figure 1: Rhombic superposition formula 



Let us denote by 9i and 82 the transforms of u>2 and u>\ via © for 
lu = u>i and 10 = 0J2, respectively. Such transformations are dehned up 
to § and, for 61 fixed, we put 9 2 = -ffli (cf. Fig.p. We have 

• the functions Q\ and 02 define the Moutard transformations of L2 
and L\, respectively, which give the same operator 

Li = —A + ui — L = L12 = —A + u\2 L 2 = —A + ?i 2 . 

• the functions ipi — ^- and (f2 — meet the equation Lip — 0. 
The general formula for u = U12 takes the form 

u = -2Alogi ^(pip 2 -P2P1) + J ((p'iP2 -pip' 2 )dz+ {pip' 2 -p\p2)dz) 

' (4) 

where u>i = pi(z) + pi(z),U2 = P2(z) + P2(z), Pi(z) and P2{z) arc 
holomorphic functions of z, and the free scalar parameter which we 
mentioned above appears as the integration constant in (QJ. 

For L — —A the Faddeev eigenfunctions © are ipo(z,X) = e Az 
and they are transformed via ^ to ipi and ip 2 - 

ipi < — Vo = e — > ip2- 

As we mentioned above these transformations are not uniquely defined 
however at least for functions lo\ and 0J2 which are polynomial in z and 
z it is possible to choose the branches of ipi and ■02 such that they would 
meet the condition |T]). Then "the cubic superposition formula" (see 
PQ and Fig. [5J implies that 

V>(z,z,A) = e Az + ^ 2 -^i) 
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Figure 2: Cubic superposition formula 



satisfies the equation Lip = 0. It seems that for smooth fast decaying 
potentials u — U12 the function ip is the Faddeev eigenfunction on the 
zero energy level. We do not state that here as a rigorous mathematical 
result however for the potentials found in [T] this is the case. 



Example. Let pi (z) = (1 



a * ~ 2' 



p 2 (z) = \{Z-bi)z 2 + ^z. 



For some appropriate constant C in 6\ we obtain 

5120|(4- i)z+ 1| 2 



(160+ |z| 2 |(4-i)z + 2| 2 ) 2 ' 

After simplifying multiplications by some constants (p± and y>2 take the 
form 

_ 2(z + I) + (4 - i)z 2 + (4 + i)z 2 
Vl ~ 160 + k| 2 |(4 — i)^ + 2| 2 ' 

_ 2(1 - i)z + 2(1 + i)z + (3 - 5i)z 2 + (3 + 5«)^ 2 
^ 2 ~ 160+ |z| 2 |(4-i)z + 2| 2 ' 

The Faddeev eigenfunction @ for £ = — A + u is as follows: 



ip(z, z, A) = e 



A 2 



1 (8i - 32)zz - 82 - (16 + 4i)z 2 - 68zz 2 
1+ A 160+ |z| 2 |(4-z)z + 2| 2 
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(32 - 8i)z + 6Sz 2 



A 2 160 + |^| 2 |(4 - + 2| s 



A 2 



as z 



The functions u, ipi, and (fi2 are smooth and real- valued, and the func- 
tion ip is also smooth in z,z. The potential u decays like l/|z| 6 as 
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\z\ — > oo, and ipi and ip2 decay like l/|z| 2 as \z\ — > oo. Hence ipi and 
if2 span a two-dimensional subspace in the kernel of L and, since u < 
and u decays sufficiently fast, there are also negative eigenvalues of L. 
We also see that "the scattering data" takes a very simple form: 

A(A,A) = --, 5(A,A) = 0. 
A 

Similar formulas for the Faddeev eigenfunctions may be easily de- 
rived for another smooth and fast decaying potential with non-trivial 
kernel which was found in [1] and corresponds to pi(z) — (i — l)z 3 + 

+ z2 +f ' P2(^) = 2iz 3 +(i + ^j) z 2 +k^z. In this case "the scat- 
tering data" are A(A, A) = -f , B(A, A) = 0. 



3 The extended Moutard transformation 
and the Novikov— Veselov equation 

In order to comply with the traditional notation related to the Novikov- 
Veselov equation we renormalize the Schrodinger operator as follows 

1 ii 

H = dd + U = -A - -. 

4 4 

In these terms the Moutard transformation of the potential and eigen- 
functions is given by the formulas 

U ^ U + 2dd logo;, 

Let us assume that ip depends also on the temporal variable t and 
satisfies the equations 

Hip = 0, d t ip = (d 3 + 5 3 + Wd + Wd)ip, dV = dU (5) 

where U = U. It is checked by straightforward computations that 

• given real-valued functions uj and V , the system |3P is invariant 
under the extended Moutard transformation [IJ 112} 

ip — > 9 = — / (ipduj — uidip)dz — (ipdui — uidip)dz+ 
u J 

[ipd 3 oj - uid 3 ip + uid 3 ip - ipd 3 uo + 2(d 2 tpduj - dipd 2 uj)~ 

2(d 2 ipduj — dipd 2 w) + 3V(ipduj — wdip) + 3V(ujdip — ipdu)]dt, 

U + 2dd\oguj, V-> \/ + 2<9 2 logw. 

(6) 
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The compatibility condition for ([5]) is the well-known Novikov- 
Veselov equation [TT] 

U t = d 3 U + d 3 U + 3d(VU) + 3d(VU) = 0, dV = dU. (7) 

It is a two-dimensional generalization of the Korteweg-de Vries equa- 
tion to which it reduces in the case U = U(x), U = V. 

We may obtain analogs of all constructions of §2 by replacing the 
Moutard transformation by its extended version (|6]) and construct 
by this mean non-trivial real-valued solutions of the Novikov-Veselov 
equation ([7]) (cf. pQ): 

• ifp\(z,t) and piiz, t) are holomorphic functions which satisfy 

dp d 3 p 
dt dz 3 ' 

then uj = pi + pi and <p = p2 + p2 satisfy JSP with U = V = 0. 
Let us consider the extended Moutard transformation of Hq = dd 
defined by w\, and let 0\ he the image of loi under this trans- 
formation. Let H = dd + U is obtained by the iteration of the 
Moutard transformation defined by 9%. Then a real-valued func- 
tion 

U(z,z,t) = 2ddlogi((piP2 -P2P1) + J {(p'iP2 ~ PiP 2 )dz+ 
(pm' ~ P\'P2)dz) + J (p1'p 2 - pip'i' + 2^ - p'(p' 2 )+ (8) 



+P1P2'" - Pi" Pi + 2(pi"p2 - PiPi"))dt ) 
satisfies the Novikov-Veselov equation 

U t = d 3 U + d 3 U + 3d(VU) + 3d(VU) =0, dV = dU. 

Example. Let pi{z) — iz 2 and p 2 = (1 + i)z + z 2 . Although these 
functions are independent on t the formula ([8} gives a non-stationary 
solution of the Novikov-Veselov equation. The Faddeev eigenfunction 
is obtained by applying iterations of (JSJ) to ipo( z >z,t) = e Xz+x *. 

These data lead to the following functions: 

TT( _ Fu(z,z,t) _ F v (z,z,t) 

U(z,z,t) = Q2j—[-^> ^')=Q27-^y 

where 

Fu(z, z, T) = Qi(2tt(iz - iz + z + z) + 2iz 4 z + Aiz 3 z - 2izz A - 4izz 3 + 
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+60iz - 60iz+ 96tzz+ 2z 4 z + z 4 + 6z 2 z 2 + 2zz 4 - 240zz- 60z + z 4 - 60z) , 
F v (z, z, t) = Qi(2At(iz - iz + z + z) + iz 4 + 4iz 3 z 2 - \2iz 2 z 3 - 6iz 2 z 2 + 
+6zzz 4 -60iz+2iz 5 +5iz 4 +60iz+48tz 2 +4z 3 z 2 +4z 3 z+12z 2 z 4 +12z 2 z 3 + 
+6zz 4 + Azz 3 - 60z - 2z 5 - 120z 2 - 60z), 

and 

Q{z,z,t) = (12-12i)t-z 3 ~(3-3i)z 2 z 2 +3iz 2 z-3zz 2 +iz 3 -(30-30i) = 

= — !— (24t - [6{x 2 + y 2 ) 2 + 8{x 3 + y 3 ) + 60]) . 
The Faddeev eigenfunctions for the potentials U(z,z,t) take the form 

ip(z,z,t,X)=e [1 + 1 ^ 

1 6(-2izz 2 - 2izz + z 2 + 2zz 2 + z 2 ) 
A Q(z,z,t) 

1 \2{iz 2 + iz - z - z 2 ) 



A 2 Q{z,z,t) 

From the latter formula it is easy to derive "the scattering data" for 
U(z,z,t): 

ip(z, z, t, A) = e Xz ^1 — + O \ T~^j J as M ~^ °°' 
i.e. we have 

A(X,\,t) = --, B(X,\,t)=0. 
A 

The function U(z, z,t) gives an example of a solution to the Novi- 
kov-Veselov equation which has smooth and fast decaying (like l/|z| 3 
as \z\ — > oo) initial data U (z, z, 0) at t = and which blows up at finite 
time (in this case it is T* = the time at which Q(z, z,t) vanishes 
at some point). The scenario of this blow up is described in pQ. 

We remark that 

1) "the scattering data" A(X, A, t) and B(X, A, t) are conserved and 
do not feel the blow-up; 

2) the real and imaginary parts of (A2(z,z,t) give solutions of the 
equation Hip = which for t = are non-singular and decay like 

O (fTp) as \ z \ ~^ 00 > i- e - they are eigenfunctions of H at the zero 
energy level. Hence the formula for ^2(2,2, t) describes a deformation 
of some functions from the discrete spectrum via the Novikov-Veselov 
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equation: these eigenfunctions stay nonsingular and quadratically in- 
tcgrable until the critical time T* when they do obtain the same sin- 
gularities as U . 
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